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ABSTRACT 

Slot timing recovery in a direct detection optical PPM communication sys- 
tem can be achieved by processing the photodetector output waveform with a 
nonlinear device whose output forms the input to a phase lock loop. The choice 
of a simple transition detector as the nonlinearity is shown to give satisfactory 
synchronization performance. The rms phase error of the recovered slot clock 
and the effect of slot timing jitter on the bit error probability were directly meas- 
ured. The experimental system consisted of an AlGaAs laser diode (X=834 nm) 
and a silicon avalanche photodiode (APD) photodetector and used Q=4 PPM 
signaling operated at a source data rate of 25 megabits/second. The mathemati- 
cal model developed to characterize system performance is shown to be in good 
agreement with the actual performance measurements. The use of the recovered 
slot clock in the receiver resulted no degradation in receiver sensitivity compared 
to a system with perfect slot timing. The system achieved a bit error probability 
of 10 -6 at received signal energies corresponding to an average of less than 60 
detected photons per information bit. 


* Work supported by the National Aeronautics and Space Administration 
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I Introduction 

Q-ary pulse position modulation (PPM) signaling is a preferred modulation 
format in free space direct detection optical communication systems. In this for- 
mat, a group of L binary source bits are transmitted as a single light pulse posi- 
tioned in one of Q=2 L possible time slots. The performance of both ideal and 
nonideal optical communication systems that use PPM signaling has been studied 
both theoretically and experimentally [l]-[6] under the condition that perfect syn- 
chronization of both the time slot boundaries and PPM word boundaries with 
those of the transmitter is available at the receiver. In practice, some means 
must be provided for the establishment of this synchronization at the receiver. 
The actual measured performance of a system that uses a phase lock loop to 
regenerate the slot clock waveform from the noisy output of a silicon avalanche 
photodiode (APD) photodetector is reported here. The system used Q=4 PPM 
signaling to transmit information at a source data rate of 25X10 6 bits/sec. The 
system was able to maintain slot clock synchronization at bit error probabilities 
less than 10 -2 , or at received signal energies that corresponded to more than 30 
detected photons per PPM light pulse. 

One of the basic problems in any communication system that uses a PPM 
signaling format is to find a means of establishing both PPM slot and word syn- 
chronization between the transmitter and receiver. The optimal receiver for a 
direct detection PPM optical communication system consists of a device which 
integrates the APD output signal over each time slot and then compares each of 
the Q integrated outputs to find the largest. Any offset or jitter in time slot 
boundaries directly affects the result of the integrations over the time slots, and, 
consequently, degrades the performance of the entire system. In practice, slot 
timing synchronization requires a regenerated clock waveform at the receiver at 
the slot clock frequency and synchronized with the clock at the transmitter. Once 
the correct slot timing is provided, PPM word synchronization can be established 
by a properly initialized modulo-Q counter driven by the regenerated slot clock. 
PPM word synchronization is maintained unless the counter miscounts or the 
recovered clock jitter is so severe that cycle-slipping occurs. Initial synchroniza- 
tion of the modulo-Q counter is usually done through recognition of a special 
data pattern transmitted at infrequent intervals. Some recent theoretical studies 
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of ways to establish PPM word synchronization are given in [7]-(9]. This work 
addresses only the issues of slot clock synchronization for receivers that use APD 
photodetectors. 

Generally speaking, all the technologies for bit timing extraction in conven- 
tional digital telecommunication systems can be used for slot timing recovery in 
optical PPM communication systems. Modification is necessary, however, to cope 
with the shot noise of photodetectors that appears with the additive Gaussian 
noise of subsequent amplifiers. The use of APD’s improves detector sensitivity 
significantly but also introduces so-called excess noise, which is nonadditive and 
follows the Conradi distribution [10][11]. Several conventional means of timing 
recovery have been studied for use in optical communications [12]- [17] and some 
of them are widely used in optical fiber telecommunication systems [18]. Free 
space optical communication systems, however, are characterized by much higher 
channel losses, typically 60-70 dB as opposed to 30-35 dB for fiber optic systems. 

This paper is organized as follows. The next section gives a brief review of 
previous studies of slot clock timing recovery for direct detection optical com- 
munication systems. Section III discusses transition detector type clock recovery 
schemes, which are widely used in conventional digital communication systems 
but have not yet been studied for optical communication systems. A theoretical 
analysis associated with this system is also developed. This method is easy to 
implement and was found to be capable of establishing satisfactory slot clock 
timing recovery. Section IV discusses the effects of slot timing jitter on system 
performance. Section V describes the experiments and measured performance of 
this slot timing recovery system. The recovered clock signal at the receiver was 
sufficient to achieve a bit error probability of 1CT 6 at a detected energy level 
corresponding to less than 60 detected photons per information bit in our 25X10 6 
bits/second laboratory system. No degradation in performance was observed in 
terms of received bit error probabilities as a function of average number of 
detected photons per information bit, when compared to the same receiver 
operated with a common transmitter/receiver slot clock (i.e. perfect receiver slot 
time synchronization). The theory developed for the design of this type of slot 
timing recovery system is shown to be appropriate as well. 
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II Review of Slot Clock Timing Recovery in Direct Detection Optical 
Communication Systems 

The response of an APD to an incident optical field can be modeled as a 
filtered compound doubly stochastic Poisson process [1, ch. 6] when the amplifier 
noise is ignored. The basic problem is to estimate slot clock waveforms from the 
times at which photons are detected from the received optical pulse train. The 
optimal procedure is to use maximum likelihood (ML) estimation. Under the 
assumptions that the APD has an unlimited frequency response, the APD output 
can be modeled as a compound Poisson random point process with intensity func- 
tion \(t) which is proportional to the intensity of the received optical field. The 
likelihood function can then be expressed as [l, ch. 3] 

PKN,; t 1 <a<t 2 } | X(t)] 

N Vs 

= exp { -J\(a)da + £ InXfo) + £ ln^rn^] } (2.1) 

ti i=l i— 1 

where pffN^; t 1 <a<t 2 } | X(t)] is the sample function density for the occurrence 
times, the uj’s, of the photon absorptions, and N ti)t2 is the total number of pho- 
ton absorptions within the interval [t 1 ,t 2 ]. The mj’s are the numbers of secondary 
electrons generated in the avalanche region of the APD in response to each pho- 
ton absorption. P(m;) is given by the Conradi distribution and does not explicitly 
contain the a v 

If the data pattern is known, i.e. X(t)==X(t-t 0 ) is known except for a time 
origin t 0 , the maximization procedure requires finding t 0 and the set of subse- 
quent pulse boundaries within [t 1 ,t 2 ] which maximize (2.1). When the data pat- 
tern is not known, the likelihood function is obtained by averaging (2.1) over all 
possible data patterns. The exact ML estimation scheme usually requires too 
many computation steps to be carried out in real time. Georghiades [7] derived 
an approximate likelihood function that reduced the amount of computation 
significantly and yet gave almost as good synchronization performance as that of 
the exact likelihood function. A digital computer is still necessary, however, to 
store the photon absorption times, evaluate the likelihood functions for different 
time slot boundary shifts, and then find the maximum through comparisons. In 
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situations where the bandwidth of the APD is limited and the data rate is high, 
the amplified APD output current cannot be modeled as a point Poisson process, 
and therefore, the ML synchronization methods., described above have to be 
modified accordingly. 

In conventional digital communication systems, a commonly used ML syn- 
chronization scheme uses an analog device to compute the derivative of the likeli- 
hood function. The result is used to derive an error signal to correct the fre- 
quency of the clock oscillator at the receiver [19]. When the data pattern is unk- 
nown, an estimated data pattern based on a previously estimated timing base 
may be used to compute the derivative of the likelihood function. This is 
referred to as data-added or decision-directed synchronization. The same method 
was studied for use in optical communication systems by Gagliardi [2]. It was 
shown that the computation becomes reasonably simple only for a few special 
pulse shapes, for example, pulses of flat top and exponential rise and fall edges. 
Rectangular shapes, which are optimal for peak power limited optical 
transmitters, such as semiconductor laser diodes, are not appropriate for this kind 
of ML estimation scheme. 

A popular variant of ML estimation for rectangular pulses is the early-late 
gate scheme [20], though it is not mathematically optimal. The early gate 
integrates the received signal over the first half of the time slot interval and the 
late gate integrates the signal over the second half of the interval, based on previ- 
ously estimated time slot boundaries. If the system is perfectly synchronized, the 
two values of the integrations will be exactly equal on average. Therefore, the 
difference of the two integrations can be used to derive an error signal to adjust 
the frequency of the slot clock oscillator. A detailed study of early-late gate syn- 
chronization in an optical communication context has been given by Gagliardi [2]. 

A phase lock loop (PLL) can also be used to regenerate the slot clock at the 
receiver. Since PPM data streams do not contain a frequency component at the 
slot clock fundamental frequency, some nonlinearity must be introduced to gen- 
erate the desired frequency component [21]. Several possible nonlinearities have 
been studied [15] [17], among them are square law, absolute value, ln(cosh) and 
fourth power. A narrowband filter can also be used to regenerate the slot clock 
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from such nonlinearly processed PPM waveforms. The use of a PLL can track 
frequency shifts of the slot clock waveform due to Doppler effects or frequency 
drift of the transmitter clock oscillator. However, the use of a narrowband filter 
results in a static phase error given by A<^ « tan _1 (2Q f Af/f 0 ), where Q f is the 
quality factor of the narrowband filter, Af is the frequency shift, and f 0 is the 
center frequency of the filter. Andreucci and Mengali [12] have studied timing 
recovery schemes that use narrowband filters for return-to-zero (RZ) data. 

Another class of timing recovery schemes, popular in conventional digital 
communication systems, are called transition detectors [20]. They generate uni- 
polar output pulses at positive or/and negative going transitions of the input 
pulses. The power spectrum of the resulting unipolar pulse sequence contains a 
frequency component at the slot clock fundamental frequency, which can then be 
tracked by a PLL or a narrowband filter. This method was used in our system 
and is demonstrated to be a simple, yet very effective means of establishing 
receiver slot clock synchronization. 

HI Slot Clock Recovery with a Transition Detector 

The shot noise of an APD is signal dependent. Unlike the noise in conven- 
tional communication systems, it is neither a stationary nor a Gaussian random 
process. If the output current of an APD is denoted by i A (t), the mean square 
value is given by [22] 

E{ii(t)> = G 2 FE{i A (t)} (3.1) 

where G is the average APD gain, F is the excess noise factor given by 
F = k eff G + (2-l/G)(l-k eff ), and k e ff is the ratio of ionization coefficients of 
holes and electrons inside the multiplication region of the APD. The average 
output current, E{i A (t)}=e\(t) is related to the received optical power as 
\(t)=? 7 P 0 (t)/hf, where e is the electron charge, tj is the quantum efficiency with 
which received photons are detected and hf is the photon energy. The APD out- 
put current is further amplified for processing in the receiver and the amplifier 
noise, which is stationary and Gaussian, has to be added to the total noise. Usu- 
ally the APD shot noise dominates when the device is operated at the optimal 
average gain for which the received bit error probability is minimum. The 
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amplified APD output signal contains mostly amplitude noise due to the multipli- 
cative nature of the APD shot noise. The transition detector, by intuition, is the 
least sensitive to amplitude noise of all the commonjy used nonlinearities that can 
be employed to introduce the desired frequency component at the slot clock fun- 
damental frequency in the received PPM waveform. Consequently, it was used 
rather than one of the other methods in our system. 


Figure 1 is a block diagram of this type of system. The input to the system 
consists of rectangular pulses of width T s corrupted by signal dependent shot 
noise and lesser amounts of additive Gaussian noise. The lowpass filter accounts 
for the limited bandwidth of the APD and the following amplifiers. The thres- 
hold crossing detector outputs a triggering pulse to the pulse shaper in response 
to a positive going transition of an input pulse. The output of the pulse shaper is 
a rectangular pulse of width T s /2. This choice of pulsewidth maximizes the 
amplitude of the fundamental frequency component at the slot clock frequency 
[20]. The pulse shaper has a ’’dead time” so that it cannot be triggered more 
than once within one time slot. This prevents false triggering by noise spikes in 
the APD output current. 


The theory of PLL tracking systems has been well developed for conven- 
tional communication systems and several studies have been done on optical com- 
munication systems as well. It has been shown that for a sufficiently narrow 
loop bandwidth, the phase error of the regenerate clock waveform has variance 
given approximately by [15] 


a 


2 

* 


2B l N(o; s ) 

Ps 


(3.2) 


Here N(o;) is the two sided power spectrum of the noise, u> s is the angular fre- 
quency of the locked PLL, P s is the power of the sinusoidal component of the 
input signal at the lock frequency, and B L is the loop bandwidth given by 

00 

B L =/ I H PLL (2rf) | 2 df, (3.3) 

0 

where Hp LL (27rf) is the transfer function of the linearized model of the PLL. 
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The phase error of the recovered slot clock obtained with a narrowband filter 
can alsq.be approximated by (3.2) with the effective loop bandwidth Bl given by 
[20] 

• 4 

oo 

B l = / I H nb ,[2ir(f - f„)] I 2 df, (3.4) 

fo 

where H nb{ (27rf) is the system function of the filter, and f 0 is the center frequency. 

Since the input to the PLL is not a stationary random process, the power 
spectrum can not be obtained by taking the Fourier transform of its autocorrela- 
tion function. The alternative is to use the definition given by [21] 

S(w) = Iim 7 ^tE{ I X T (w) I =} (3.5) 

T— co 21 

where Xx(cu) is the Fourier transform of a sample of the signal waveform from -T 
to T given by 

T 

X T (cu) = / x(t)e^ wt dt (3.6) 

-T 

As in Chen’s approach in [15], the input signal to the PLL is defined as 

x(t) = £ p (t-nQ T s -c n T s -r n ) (3.7) 

n 

where p(t) represents the pulse shape, T s is the duration of a PPM time slot, c n is 
an integer from 0 to (Q-l) used to indicate the random PPM pulse position, and 
T n represents random jitter in threshold crossing times. If all transmitted data 
are equally likely, the c n ' s form a set of independent, identically distributed 
(i.i.d) random variables. The r a ' s are also considered to be i.i.d. random vari- 
ables. Appendix A shows that under these assumptions, the power spectrum (3.5) 
may be expressed as 

SM = | pm 1 2 ( JL-(i- 1 CM | 2 1 m,m | 2 ) + 

+ I M,M | 2 g «M2>*/T S )|, (3.8) 

(QTs)^ k=-oo 

where P(cu) is the Fourier transform of the pulse shape p(t), C(c o) and M^cu) are 
the characteristic functions of the c n and r n respectively. These are given by 
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C(oj) = E{e iwTsC "} = -J-^ 1 e jwTsn (3.9) 

Q n=0 

and 

4 4 

M^w) = E{e jwrn }. (3.10) 

The first term in (3.8) corresponds to noise and the second term corresponds 
to the discrete frequency components at the slot clock fundamental frequency and 
its harmonics. Therefore, the noise power spectrum may be written as 

N(w) = -±- I PM 1 2 ( 1 - I CM I 2 | MrM I 2 ). (3.11) 

The total power at the slot clock fundamental frequency is given by 

P. = — I MrM) I 2 I PM) I 2 - (3-12) 

If the random jitter represented by the r n ’s can be considered as zero mean 
Gaussian random variables, then 



M t (w) = e~ 2 ' (3.13) 

where <7 r 2 represents the variance of the threshold crossing times which will be 
determined next. 

The input signal to the threshold crossing detector consists of the APD out- 
put current, i A (t), and amplifier thermal noise, n t (t). The threshold crossing time, 
t c , is the solution to the equation i A (t c ) + n t (t c ) = 1^, where I th represents the 
threshold. The APD output current, i A (t), can be expressed as the sum of the 
signal and shot noise, as 

i A (t) = s(t) + n s (t) (3.14) 

where s(t) = E{i A (t)} and n s (t) = i A (t) - E{i A (t)}. At high input optical signal 
levels, the threshold crossing time, t c , varies within a small region about its mean 
value T c . Consequently, it is appropriate to approximate s(t) by the first two 
terms of its Taylor expansion about T c , and write 

Ith ~ S (T C ) + s' (T c )(t c -T c ) + n s (t c ) + n t (t c ). (3.15) 
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Since n s (t) and n t (t) are independent and s(T c )=I th , it follows that 

[s' (T c )l V = E{n, 2 (t c )} + E{n t 2 (t c )}. 

The variance of the amplifier noise is given by 

4kT r B 


EMtJ} = 


R 


(3.16) 


(3.17) 


where k is Boltzmann’s constant, T r is the amplifier equivalent noise temperature, 
B is the bandwidth of the lowpass filter, and R is the load resistance seen by the 
APD. 

The variance of the shot noise, E{n s 2 (t c )}, is equal to the variance of the 
APD output current, which is a function of t c . Conditioned on t c , the APD out- 
put current is a filtered compound Poisson process with variance given by [l, ch. 

4] 

tc 

E{n s 2 (t c ) | tj = e 2 G 2 F f X(t)h 2 (t c -t)dt. (3.18) 

—00 

where h(t) is the impulse response of the lowpass filter. \(t) is the photon count- 
ing intensity function given by 


X(t) = 



t>t 0 
t <t Q 


(3.19) 


where t 0 is the pulse transition time prior to t c . If the laser has a high on-off 
extinction ratio and the background radiation appears to be small, i.e. X 1 »X 0 , 

te 

E(v(t c ) | tj « e-FG 2 \i/h 2 (t c -t)dt 

to 

T c -to t c -to 

= e 2 FG 2 X 1 [ / h 2 (u)du + / h 2 (u)du]. (3.20) 

o Tc-to 


Since t c -T c is very small, the integrand of the second term in (3.20) can be 

approximated by h 2 (u)^h 2 (T c -t 0 )+2h(T c -t 0 )-^-h(T c -t 0 )[u - ( T c - t 0 )). Averaging 

dt 

both sides of (3.20) over t c , 

T«rt° 

E{n s 2 (t c )> « e^G^X, [ / h 2 (u)d<i + h(T„-g-ih(T c -t 0 K 2 ]. (3.21) 
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Next, according to [1, ch. 4], 

t 

s(t) = E{i A (t)} = eG f X(u)h(t-u)du (3.22) 

-00 


then 


s' (t) ~ eGXih(t-t 0 ) 


(3.23) 


where X!»X 0 is assumed. Substituting (3.17), (3.21) and (3.23) into (3.16), the 
variance <7 r 2 can be solved as 


Tc-to 

FG 2 Xj / h 2 (u)du + 
o 


4kTrB 

Re 2 


[GX.hfT.-tof - FG 2 X,h(T t -t 0 ) Ah(T e -t 0 ) 


(3.24) 


Equation (3.24) may be evaluated by assuming that the lowpass filter is ideal, i.e. 

/ e " j25rftd I f I <B 

H(j27rf) = \ Q 0therwi5e . (3-25) 


and 


h(t) = 2B 


sin2;rB(t-t d ) 

27rB(t-t d ) 


(3.26) 


where t d is the time delay of the filter. 

There is an optimal point for threshold crossing time T c , at which (3.24) is 
minimum. In practice, it corresponds to an optimal threshold level applied to the 
threshold crossing detector. Figure 2 is a plot of the normalized rms threshold 
crossing time jitter, a T /T s , as a function of normalized mean threshold crossing 
time, ATcB, where AT C = T c -(t 0 +t d ). In other words, AT C is denoted as the 
time difference between the mean threshold crossing time, T c , and the midpoint 
of the rising edge of the input pulse, t 0 +t d . The curve in figure 2 is nearly flat 
near BAT C = -0.1. In fact, a r j T s increases by only 3% when BAT C changes from 
its optimal value ( «-0.08 ) to zero. Consequently, the value BAT c =0 can be 
used in (3.24) to a good approximation. This corresponds to a threshold level at 
one half the peak amplitude of the APD output pulses. As shown in Appendix B, 
when BAT c =0, (3.24) reduces to 
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a 


2 

r 


FG 2 Xj + 


4kT r 

Re 2 


4BG 2 Xj 2 


(3.27) 


Finally, substituting (3.27), (3.13), (3.12), (3.11) and (3.9) into (3.2), the vari- 
ance of the phase error in the recovered slot clock is given by 


= QTsB L *{exp [ 


*2 FG 2 XjT 3 + 4kT r T s /Re 2 


BT C 


G 2 (X 1 T s ) 2 


]-l}. 


(3.28) 


The phase error variance, a decreases as the bandwidth of the lowpass 
filter, B, increases. However, the effect of false triggering by noise spikes is not 
included in (3.28). This may contribute significantly to the overall phase error 
when the bandwidth of the lowpass filter becomes too large. In practice, B is also 
limited by the frequency responses of the APD and the amplifiers. In our experi- 
ment, B=2/T s was used. 


IV Effects of Slot Clock Jitter on Probability of Bit Error 

Jitter in recovered slot clock timing at the receiver results in the appearance 
of part of the received PPM pulse energy in an adjacent PPM slot, and hence 
acts as a source of background noise. If the integrators are ideal, the amount of 
signal energy which spills over into the adjacent slot is proportional, on average, 
to the amount of offset in the slot boundaries. If the timing offset of the 
recovered slot clock is r, the fractional offset is given by e = r/T s , and the 
effective average number of signal photons becomes 

X,' T s = (l-e^T,. (4.1) 

The effective average number of noise photons in the adjacent slot becomes 

\q T s = X 0 T S + eXjT s . (4-2) 

The average number of noise photons in the remaining Q-2 time slots is still 
given by X 0 T S . The probability of a PPM word error (PWE) is now given by 

OO X X 

PWE(c) = 1 - / P(x | X,' T,)/P(y | X 0 ' T,)[/P(z | XoTJdzJ^dydx (4.3) 

0 0 0 

where P(u | XT S ) is the probability density function for the APD output given 
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that the average number of absorbed photons was X/ T s , X 0 ' T s , and X 0 T S 
respectively. 

Equation (4.3) usually requires excessive computation when evaluated 
numerically due to the complicated distribution of the APD output. The alterna- 
tive is to use the union bound as shown below, 

00 00 

PWE(<) « fP(x I X,' T s )/P(y I X 0 ' T s )dydx 

0 x 

00 00 

+ (Q-2)/P(x I X,' T s )/P(y I X„T s )dydx. (4.4) 

0 x 


The union bound gives a good approximation when the input signal to noise ratio 
is high and the timing offset is small, i.e. X 1 /X 0 »l and e«l. A detailed 
description of an efficient numerical evaluation of (4.4) is given in [6]. The 
overall probability of a PPM word error is obtained by averaging (4.4) over e. If 
the timing jitter is assumed to have Gaussian distribution with zero mean and 
standard deviation given as o t — it. follows that 


CO 

PWE = / PWE(e) 
-00 


1 

\/27rcr £ 



(4.5) 


The probability of bit error (PBE) is given by PBE 


_ Q 

2 ( 0 - 1 ) 


•PWE [23]. 


In reality, none of the high speed integration circuits implemented in 
hardware is ideal and the amount of received signal energy in adjacent slots can 
not be determined simply by using (4.1) and (4.2). As an example, the integra- 
tions over the PPM time slots in our experiment were realized by sampling the 
output of the tapped delay line matched filter shown in figure 3. Ideally, the 
impulse response of the matched filter should be a rectangle of width T s , so that 
the output in response to a rectangular input pulse is a symmetric triangle of 
base 2T S . A sample at the top of the triangle is proportional to the integration of 
the energy contained in the input signal within the previous time interval of 
length T s . There is no intersymbol interference because the triangle waveform 
vanishes completely at the next sampling time. However, due to the limited 
bandwidth of the system, all three corners of the triangle become rounded. If the 
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input rectangular pulse has zero rise time and fall time, and the lowpass filter has 
a "brick wall” like frequency response with cutoff frequency 2/T s , the output of 
the matched filter shown in figure 3 can be obtained through Fourier analysis. 
Figure 4 is a plot of such a distorted triangle along with the associated perfect 
triangle (dashed line). In practice, input rectangular pulses have finite rise and 
fall times, and no lowpass filters have true "brick wall” frequency response. It 
may be impossible to derive the actual distorted triangle shape, and consequently 
determine the exact amount of received signal energy spilled over into an adja- 
cent slot due to timing offset. When the triangle output by the matched filter 
becomes rounded, the values of samples at these turning points of the triangle 
become less sensitive to the sampling times. Therefore, jitter in the recovered slot 
clock has less effect on the system performance (PWE and PBE) than that when 
a perfect integrator is vised. Nevertheless, (4.5) can be used as an upper bound 
for the effect of jitter in the recovered slot clock on system performance. 

V Measurements of System Performance 

A prototype optical communication system was built that consists of a GaA- 
lAs diode laser (Sharp LT024MD, \=834nm) as the transmitter and a silicon 
APD (RCA C30902S) as the photodetector. The system used Q=4 PPM signaling 
at a source data rate of 25 Mbits/second. A detailed description of the system is 
presented in [5] and [6]. The slot clock recovery system shown in figure 1 was 
constructed. It used a high speed comparator (Motorola MC1650) as the threshold 
crossing detector and a monostable multivibrator (Motorola MC10198) as the 
pulse shaper. The inherent recovery time (dead time) of the multivibrator 
prevented it from triggering more than once within one slot period. The PLL 
consisted of a double balanced analog mixer (Motorola MC12002), an active 
second order loop filter, and a voltage controlled crystal oscillator (VCXO). The 
effective loop bandwidth was designed and experimentally verified to be about 
B l =1 KHz. 

Timing errors were measured with the use of an rms phase error detector 
shown in figure 5. A double-balanced frequency mixer (ASM-15) was used as the 
phase error detector. Under the conditions of small phase error, the output 
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voltage can be approximated as V pd =K pd A0, where K pd is the detector gain, and 
Ad is the phase difference between the two input signals. The amplified output 
of the phase detector was input to an rms to DC,. converter (Burr-Brown 4341) 
which had a frequency response from DC to 450 kHz. The output of the phase 
detector often contained a DC component because a static phase shift between 
the two input signals was difficult to balance out completely. Since the phase 
jitter was often very small compared to static phase shifts, a coupling capacitor 
was inserted between the phase detector output and rms to DC converter input. 
The coupling capacitor and the input impedance of the rms to DC converter 
formed a high pass RC filter. A large coupling capacitor was chosen so that the 
3dB cutoff frequency was below 0.1 Hz. Two identical slot clock waveforms with 
a series of known fixed phase shifts were used to calibrate the circuit. It was 
found that the circuit had a nearly constant gain of 7.62 V/radian for input 
phase shifts up to ±tt/ 4 radians. 

Other sources of phase errors that appear in the phase error measurements 
include noise of circuit components, such as the amplifiers, and the small phase 
jitter inherent in the VCXO. These phase errors were not included in our model 
used to compute the phase error, and therefore, should not be included in the 
measurements. The rms value of these phase errors was determined by feeding 
the jitter free clock signal into the PLL and then measuring the phase error of 
the VCXO output against the input jitter free clock. Since the phase jitter of the 
recovered slot clock is statistically independent of the VCXO noise and other cir- 
cuit noise, the following relationship holds, 

^total = &<p "h ^<pn 

where ^totai is total measured rms phase error, is the rms phase error of 
the recovered clock described by (3.28), and cr^ n accounts for the phase error due 
to the noises of the VCXO and the circuit. The jitter of the recovered slot clock 
is then given by 


a 4> 



^(j) total [ 1 


1 , ^ 0 ' 

~ total 


(5.2) 


Measurements revealed that <7^ n was about 7.6 X10 -4 radians. The smallest value 
of c^total measured was about 5X10 -3 , so that the approximation a ^ ~ c^total 
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was used throughout the measurements. 


The optimal level for the threshold crossing detector was set experimentally 
such that the observed phase error achieved a minimum. As predicted in section 
HI, the optimal threshold was found to be near one half of the average peak 
amplitude of the APD output pulses, and small deviations from the optimal level 
resulted in negligible increases in the measured rms phase error. Figure 6 is a 
plot of the normalized rms phase error of the recovered slot clock, a^/ 2 k, as a 
function of the average number of detected photons per information bit given by 


photons/bit = 


VTs_ 

log 2 Q 


*?P a yQT s 

hflog 2 Q 


(5.3) 


Here log 2 Q is the number of information bits contained in a PPM word, tj is the 
quantum efficiency of the APD, taken to be 77%, and P av is the average optical 
power incident on the APD. P^v was directly measured by placing an optical 
power meter at the position of the APD. The measured phase errors of the 
recovered slot clock using a crystal bandpass filter of effective loop bandwidth 
B l =9.5 KHz (full bandwidth 19KHz) is also plotted in figure 6. The solid and 
dashed curves were obtained from the theoretical analysis. They agreed very well 
with the measurement data. 

Finally, the recovered slot clock was used in the receiver and the system bit 
error probability was measured as a function of the average number of detected 
photons per information bit. The data source consisted of pseudo random binary 
sequences, 1023 bits long. The regenerated binary sequence at the receiver was 
compared with the properly delayed source binary sequence with the use of an 
exclusive OR gate. The output of the exclusive OR gate was sampled at the 
source data rate to form a series of short pulses each of which corresponded to a 
bit error. A counter was then used to counts the bit errors, and the probability 
of bit error was obtained by dividing the bit error rate by the source bit rate. 

Figure 7 shows the results obtained with the slot clock regenerated by the 
PLL (B l =1000Hz) along with the results obtained with the jitter free slot clock 
(the clock common to the transmitter) at the receiver. The solid curve 
corresponds to the numerical computations described in [6] that assumed perfect 
slot timing, and the dotted curve represents the numerical results for the upper 


A 
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bound given in (4.5) with determined by (3.28). It is shown that the recovered 
clock worked so well that no penalty was observed in terms of detected photons 
per bit at a fixed bit error probability. No cycle slipping and losses of lock were 
observed until the received optical power dropped below 15 photons per bit (30 
photons per pulse) which corresponded to a bit error probability above 0.01. Fig- 
ure 8 is similar to figure 7 except that the bandpass filter (B L =9.5KHz) was used 
in place of the PLL. The dotted curve is again the upper bound given by (4.5) 
and (3.28). Cycle slipping in this case occurred only when the detected optical 
power dropped below 30 photons per bit (60 photons per pulse). 

Our results for the rms phase errors of the recovered slot clock are compar- 
able with the theoretical results obtained by Chen [15] and Ling [16], which used 
different methods but did not consider the effect of the APD excess noise and 
amplifier noise. For example, under the conditions of B L T S = 2.0X10 -5 , 60 signal 
photons per information bit, and no background radiation, the result in [15] 
corresponded to cr^/2^ =3.5 X1CT 4 using a square law device followed by a PLL. 
The result in [16] with an early-late gate scheme was given by 0^/2^ =8.2X10 -4 . 
Our result under the same conditions was <7^/ 27 T =9.8X10~ 4 , which included the 
effects of the APD shot noise and amplifier noise. 

VI Conclusion 

We have developed, both in theory and experiment, a transition detector 
PPM slot clock recovery system. It is easy to implement and gives excellent per- 
formance. An upper bound for the probability of bit error with the use of the 
recovered slot clock in the receiver was derived. The measured performance was 
shown to be below the upper bound. The limited bandwidth of the matched 
filter made system performance less sensitive to small amounts of jitter in the 
times at which the output waveform of the matched filter were sampled. 



Appendix A 

Power Spectrum of PPM Signal 


- 18 - 


Substitute (3.7) into (3.6) and let T=NQT S , , 

X t H = £ P(w)e" jw(nQT>+CnT ‘ +ril) . (A.l) 

n=-N 


Since the c n ’s and the r n ’s are two independent sets of i.i.d. random variables, 

E{ | X t (cj) | 2 ) = | P(w) | 2 f) e j ^ T ^ m " n ^{e jwT ^ Cm_Cn) }E{e jw(rm ' rn) } 

m,n=-N 

= I P(w) I 2 [ (2N+1) + 

+ | E{e jwTsCn } | 2 | E{e JWTn } | 2 £ e jwQTs(m " n) ] 

m^n 

= I P(o>) I 2 [ (2N+1) (1 - I CM I 2 I MM) I 2 ) 

+ | C(w) | 2 | MM) I 2 I E I 2 I- (A.2) 

n=-N 

The power spectrum (3.5) can be obtained by taking the the limit of (A.2) as 
T— >oo, therefore, 

S(w) = | PM | 2 | 1- | CM I 2 I MM) I 2 ) + 

+ I CM I 2 I MM) I 2 Jim ± | £ | % (A. 3) 

T— oo 2 I n== _N 


Since e _jXn = 2ir <5(x- 27rk), and <5(ax)=- — — <5(x), then, 


n=-oo 


k=-oo 


N 


1 S N ^ 1 2 - X 1 ” wl ^ 2 ™/« T >)l 2 


1 , 27T 


00 


n — -oo 


4n ' E Jim ^r^im/QT,). 


(QT s ) 2 „iJooT-o= 2T 


(A.4) 


Using the relationship that 5(u>)=-^— lim 2T S — ■ - , each term in the sum of 
" 27TT-00 wT 

(A.5) reduces to 

1 _ 11 sin(a;-27m/QT s )T 

l ° Lw '» r 2T m 2 wq T ,)t 
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= --^<5(w-27rn/QT s ). (A.5) 


According to (3.9), 



1 CM | = A 

1 _ e jwT)Q . 1 . sin(wT s Q/2) 

(A.6) 


1 i.eJwT. 1 Q 1 sin(o-T s /2) 1 ’ 

therefore, 

| C(2jrn/QT S ) | 

=1 if n=kQ for any integer 

k, and 


| C(27rn/QT S ) | =0 otherwise. Substituting (A.4), (A.5), and (A.6) into (A.3), the 
power spectrum of PPM signal is given by 

S(w) = I p (“) I 2 [ I CM | 2 1 M^) I 2 ) + 

+ I MM I 2 E %-- Irk /T.)l. (A.7) 

(QT s f k=-oo 


Appendix B 

Variance of Threshold Crossing Times at BAT c =0 

The impulse response of the low pass filter (3.26) and its derivative at time 
t=T c -t 0 can be rewritten by substituting the new variable AT c =T c -(t 0 +t d ), as, 

sin27rBAT„ 


h(T c -t 0 ) = 2B- 


27tBAT„ 


rl 27rBAT e cos27rBAT c - sin27rBAT c 

-2-h(T c -t 0 ) = 2B — r 

dt ^ c 2 ttBAT 2 


(B.l) 

(B.2) 


Next, 


T c -to 


Tj-to-td 


AT, 


J h 2 (u)du = / h 2 (v+t d )dv « J (2B 


-t d 


-00 


sin27rBv 

27rBv 


) 2 dv 


2xBAT e 

= B+® / 


7T 


(B.3) 


kj m 

where the relationship J ( Sinx -) 2 dx = has been used. When BAT C — ►(), equa- 


-oo 


tion (B.l), (B.2), and (B.3) reduce to h(T c -t 0 ) = 2B, -^-h(T c -t 0 ) = 0, and 
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TV to 

J h 2 (u)du = B. Substituting these into (3.24), the variance of threshold cross- 
o 

ing time at BAT c =0 becomes 


a 


2 

r 


FG 2 Xj + 


4kT r 

Re 2 


4BG 2 X 1 2 


(B.4) 
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Figure Captions 

Fig. 1. Block diagram of a slot clock recovery system with a transition detector. 

Fig. 2. Normalized rms threshold crossing time, cr r /T s , vs. normalized mean thres- 
hold crossing time, BAT C . 

Fig. 3. Block diagram of the matched filter. 

Fig. 4. The output of the matched filter shown in Fig. 3 in response to a rec- 
tangular input pulse. The dotted line corresponds to the output of an ideal 
matched filter. 

Fig. 5. Block diagram of the rms phase error detector. 

Fig. 6. Normalized rms phase error of the recovered slot clock, o ^/ 2 - k , vs. average 
number of photons per bit. The curves are obtained by numerical evaluations of 
(3.28). The squares and triangles represent the experimental measurements. 

Fig. 7. Probability of bit error vs. average number of photons per bit when the 
slot clock regenerated by a PLL (B l =1 KHz) was used. The background noise 
intensity was n b =0. 0122/slot. The solid curve is obtained as in [6] assuming per- 
fect slot clock. The dashed curve represent the upper bound given by (4.5). The 
small crosses and triangles represent the experimental data. 

Fig. 8. Same as Fig. 7, but with the slot clock regenerated by a narrowband filter 
(B l = 9.5 KHz). 
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